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NOMENCLATURE 


B  average  thickness  of  the  surface  frazil  layer 

B(-  critical  surface  layer  thickness  for  frazil  structure  development 

Bniax  thickness  of  surface  layer  with  all  frazil  collected  on  the  surface 

C  local  vertical  frazil  concentration  profile  below  the  surface  layer 

C  mean  volumetric  frazil  concentration  in  the  water  column  below  the  surface  layer 

c  average  or  global  vertical  concentration  profile  below  the  surface  layer 

C(i  drag  coefficient 

Q  frazil  concentration  of  the  surface  layer 

C*  dimensionless  Chezy  conveyance  coefficient  of  a  river  reach 

/  frequency  of  the  surface  renewal 

ft  dimensionless  average  surface  age 

ft*  dimensionless  surface  renewal  frequency 

fi  dimensionless  time 

F  probability  distribution  function  for  averaging  local  unit  areas  over  the  water  surface 
g  acceleration  due  to  gravity 

(l-Pi/P)^ 

H  flow  depth 

k  von  Karman  constant 

L  latent  heat  of  fusion 

M  mass  of  frazil  in  the  water  column  under  a  unit  surface  area 
N  number  of  crystals  per  m^ 

Q  heat  loss  from  the  water  surface 
r  mean  radius  of  the  frazil 

R  Reynolds  number  of  rising  frazil  crystals 
t  time 

i  average  surface  age 

/*  time  for  a  crystal  to  rise  from  the  bottom  of  the  flow  to  the  surface 
time  of  initiation  of  equilibrium  growth  of  frazil 
u*  shear  velocity 

U  minimum  flow  velocity  to  maintain  well-mixed  frazil  in  a  river 
V  terminal  rise  velocity  of  the  frazil 

w  frazil  crystal  thickness 

y  distance  from  the  bottom  of  the  flow 

a  diameter/thickness  ratio  of  frazil  crystals 
P  mass  fraction  of  the  frazil  in  circulation 

p  density  of  water 

Pi  density  of  ice 

T  time  from 
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A  Model  for  Vertical  Frazil  Distribution 


CHYR  PYNG  LIOU  AND  MICHAEL  G.  FERRICK 


INTRODUCTION 

Frazil  ice  consists  of  individual  disk-shaped  crystals  formed  in  turbulent  water.  It  is  initiated  by 
seeding  at  the  air/water  interface,  followed  by  secondary  nucleation  (Osterkamp  1978.  Osterkamp  et 
al.  1983,  Daly  1984,  Ashton  1986).  Wide  ranges  of  frazil  disk  diameter  and  thickness  have  been 
reported  (Schaefer  1950,  Arakawa  1954,  Kivisild  1970,  Daly  and  Colbeck  1986).  Due  to  the  higher 
growth  rate  in  the  basal  plane  of  an  ice  crystal  (Hobbs  1974).  young  frazil  crystals  are  smaller  and 
thicker  than  older  ones.  Frazil  crystal  concentrations  in  natural  rivers  range  from  10  thousand  to  1 
million  crystals  per  cubic  meter  (Schaefer  1950,  Osterkamp  and  Gosink  1982).  If  we  assume  a 
diameter/thickness  ratio  of  30,  the  volume  concentration  corresponding  to  1  million  3-mm-diam. 
crystals  per  cubic  meter  is  0.0007;  thus  the  frazil-water  mixture  is  very  dilute.  Frazil  crystals  that  are 
actively  growing  in  slightly  supercooled  water  can  attach  to  any  solid  surface  upon  contact,  and  can 
lead  to  ice  dams,  anchor  ice,  and  icing  of  water  intake  structures  in  rivers.  Frazil  is  also  believed  to  be 
an  important  mode  of  ice  growth  in  the  marginal  ice  zones  of  the  Arctic  and  the  Antarctic.  Osterkamp 
(1978),  Martin  (1981),  Daly  (1984)  and  Ashton  (1986)  review  frazil  formation  and  growth. 

The  evolution  of  frazil  into  larger  surface  ice  forms  such  as  floes,  pancakes,  and  eventually  floes  is 
not  completely  understood.  Michel  (1978)  proposed  that  frazil  crystals  agglomerate  to  form  porous 
floes  within  the  flow  during  the  active  growth  stage,  and  then  concentrate  at  the  water  surface  as  frazil 
slush.  If  the  slush  is  dense  and  stays  at  the  surface  long  enough,  ice  floes  in  the  shape  of  pancakes 
develop.  Larger  surface  ice  forms  may  evolve  from  these  pancakes.  However,  since  the  frazil-water 
mixture  is  very  dilute,  the  chance  for  crystals  tocollide  and  form  floes  in  the  flow  interior  is  low  relative 
to  that  at  the  surface,  and  has  not  been  well  documented. 

Because  they  are  buoyant,  frazil  crystals  tend  to  concentrate  near  the  water  surface,  increasing  the 
opportunities  for  collisions  between  crystals.  The  crystals  at  the  water  surface  in  direct  contact  with 
subfreezing  air  have  a  greater  potential  for  heat  loss  and  freeze  bonding  than  those  at  depth.  Therefore, 
structured  frazil  accumulations  are  more  likely  to  develop  at  the  water  surface  than  in  the  interior  of 
the  flow.  Martin  and  Kauffman  (1981)  found  from  their  saltwater  wave  tank  experiment  that  frazil 
crystals  floated  to  the  water  surface  and  were  herded  by  the  waves  to  one  end  of  the  tank.  The  individual 
crystals  rapidly  sintered  together  to  form  large  irregular  clumps.  Osterkamp  and  Gosink  (1982) 
suggested  that  at  high  velocities,  turbulence  may  prevent  flocculation  and  ice  pan  formation,  and 
instead,  frazil  slush  may  develop  at  the  surface.  If  the  turbulence  is  strong  enough,  frazil  may  be  kept 
in  suspension  and  the  water  surface  may  be  essentially  free  of  ice.  Hirayama  (1986)  reported  field 
observations  of  both  suspended  and  layered  frazil.  Tsang  (1988)  proposed  a  model  for  frazil 
distribution  based  on  energy  considerations,  and  compared  the  results  with  limited  field  data. 

In  this  report  weconsider  the  processes  in  a  turbulent  flow  that  determine  if  frazil  crystals  will  remain 
in  suspension  or  if  more  structured  ice  will  form  at  the  water  surface.  We  propose  a  model  to  enhance 
our  understanding  of  the  physical  proces.ses  involved  in  frazil  evolution,  to  systematically  assess 
observations  of  mixed  and  layered  flows,  and  to  provide  a  logical  structure  to  guide  further  work. 


FRAZIL  GROWTH 


Frazil  grows  through  the  removal  of  the  heat  of  fusion  by  turbulence  and  the  eventual  heat  transfer 
to  the  air  above  the  water  surface.  Carstens  ( 1 966)  and  Osterkamp  (1978)  present  an  idealized  water 
temperature  history  during  periods  of  frazil  production  and  growth.  When  the  air  temperature  is 
subfreezing,  the  water  temperature  can  drop  below  the  freezing  point.  Frazil  nucleation  begins  soon 
after,  and  the  rate  of  temperature  decrease  slows  as  latent  heat  is  released  by  the  growing  crystals.  The 
water  temperature  then  reaches  a  minimum  and  begins  to  recover  as  the  rate  of  ice  growth  increases. 
Eventually  an  equilibrium  is  reached  at  time  when  the  rate  of  heat  released  by  the  crystals  equals  the 
rate  of  heat  loss  to  the  air.  The  dynamics  oi  initial  frazil  formation  and  growth  are  not  well  understood. 
The  modes  of  ice  growth  that  can  occur  are  crystal  number  growth  and  crystal  size  growth.  The  initial 
period  of  nonequilibrium  growth  prior  to  is  believed  to  be  brief  (minutes).  After  t^,  a  slight  amount 
of  supercooling  in  the  water  is  maintained,  which  is  responsible  for  the  continued  growth  of  frazil 
crystals.  Detailed  models  of  nonequilibrium  frazil  growth  can  be  found  in  Daly  (1984)  and  Omstedt 
(1985). 

During  the  equilibrium  growth  period  after  the  time  rate  of  increase  of  frazil  mass  M  in  the  water 
column  under  a  unit  surface  area  can  be  related  to  the  rate  of  heat  loss  Q  from  the  surface  by  the  principle 
of  energy  conservation; 


L^=Q(x) 

dz 


(1) 


where  L  is  the  latent  heat  of  fusion  and  x  denotes  time  since  Equation  1  is  integrated  to  obtain  the 
total  mass  of  frazil  in  the  water  column  at  x; 

M(x)=lf  Q[t)dt  (2) 

^  Jo 

The  integration  constant  can  be  determined  using  Af  =  A/f  at  x  =  0.  The  ice  mass  increases  with  time 
as  long  as  the  surface  heat  loss  continues. 


TURBULENCE  AND  FRAZIL  TRANSPORT 
NEAR  THE  WATER  SURFACE 

Large  surface  features  such  as  boils  and  eddies,  attributable  to  turbulence,  are  often  visible.  Komori 
et  al .  ( 1 982)  conducted  a  laboratory  study  on  turbulence  and  transport  mechanisms  at  the  water  surface 
in  a  shallow  open-channel  flow.  They  found  that  1)  the  intensity  of  vertical  velocity  fluctuations  is 
reduced  near  the  free  surface,  which  increases  the  intensity  of  velocity  fluctuations  in  the  streamwise 
and  lateral  directions,  2)  large  energy-containing  eddies  well  up  toward  the  surface,  spread  horizon¬ 
tally,  then  plunge  back  into  the  flow  interior,  and  3)  the  time  fraction  when  these  eddies  are  at  the 
surface  is  only  about  10%,  but  they  contribute  more  than  60%  of  the  mixing  near  the  surface.  Similar 
turbulence  structures  with  60%  of  the  Reynolds  stress  produced  by  large  eddy  motions  occupying  only 
10%  of  the  time  were  observed  by  Seitz  (1973)  in  an  estuary  and  by  Gordon  and  Dohne  (1973)  in  a 
river.  Additional  information  on  this  subject  can  be  found  in  Brown  (1970).  Davis  and  Lozano  ( 1 979). 
Brutsaert  and  Jirka  (1984),  Shimko  (1985)  and  Sonin  et  al.  (1986). 

We  have  also  observed  large  eddy  motion  during  experiments  in  the  36-m-long  and  1.2-m-wide 
tilting  flume  at  the  Cold  Regions  Research  and  Engineering  Laboratory.  A  flow  rate  of 0.057  m  Vs  was 
used  with  both  a  bare  aluminum  channel  floor  and  a  roughened  channel  floor  that  increased  the 
Manning’s  roughness  of  the  flume  by  a  factor  of  10.  Talc  powder  was  spread  over  the  water  surface 
near  the  upstream  end  of  the  flume,  and  this  patch  of  the  surface  was  videotaped  with  a  camera  mounted 
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on  an  instrument  carriage  that  moved  with  the  flow.  The  upwelling  of  interior  water  to  the  surface  was 
vivid  as  the  eddy  motions  temporarily  cleared  the  talc  at  various  locations  in  a  seemingly  random 
pai  .cm.  The  sizes  of  the  renewed  areas  were  of  the  order  of  the  How  depth.  This  renewal  behavior  was 
.  s  1  for  all  flow  conditions  except  for  the  lowest  velocity  (0. 1  m/s)  case  in  the  smooth  channel. 

Wc  piopose  that  the  surface  renewal  eddies  act  as  the  primary  mechanism  whereby  frazil  crystals 
at  the  water  fa  :e  are  entrained  and  redistributed  over  the  depth  of  the  turbulent  flow.  A  large  eddy 
; 1  the  wau  urface  that  is  unifonnly  covered  with  frazil  crystals  and  momentarily  clears  an  area. 
.As  iiie  eddy  plunges  back  into  the  flow  interior,  it  entrains  some  of  the  crystals  that  occupied  the 
temporarily  cleared  area.  We  assume  that  the  entrained  crystals  and  the  crystals  in  the  flow  interior  are 
unifonnly  redistributed  over  the  flow  depth  by  the  eddy  mixing.  Afterward,  the  net  upward  motion 
resumes.  The  large  eddies  also  promote  horizontal  surface  currents  as  they  spread.  The  opposing 
surface  currents  between  closely  spaced  eddies  herd  the  frazil  crystals  into  the  narrow  region  between 
them.  The  concentrated  crystals  are  exposed  to  subfreezing  air  and  may  agglomerate  to  form  floes. 

The  frequency/of  the  surface  renewal  affects  the  time  available  for  the  frazil  at  the  water  surface 
to  become  more  structured.  We  assume /to  be  uniform  over  the  water  surface  and  constant  over  time. 
Methods  of  estimating/can  be  found  in  the  literature  on  gas  absorption  over  turbulent  liquid  surfaces. 
O’Connor  and  Dobbins  (1958)  present  an  argument,  based  on  mixing  length  theory,  that/equals  the 
vertical  velocity  gradient  at  the  surface.  By  assuming  a  logarithmic  velocity  distribution,  they  showed 
that 


where  k  =  the  von  Karman  constant 
H  =  the  flow  depth 
u*  =  the  shear  velocity. 

Komori  et  al.  (1982)  measured/to  be  0.82  Hz,  which  expressed  in  terms  of  shear  velocity  and  flow 
depth  replaces  X/k  in  eq  3  by  3.73.  We  will  use  eq  3  with  the  experimentally  determined  constant  to 
estimate  /. 


FRAZIL  TRANSPORT  IN  THE  FLOW  INTERIOR 


Between  renewals,  the  rise  velocity  due  to  the  buoyancy  of  frazil  in  the  interior  depends  on  size, 
shape  and  orientation  of  the  crystals,  and  on  the  turbulence  associated  with  eddies  smaller  than  the 
energy-containing  ones.  Since  frazil  crystals  have  a  range  of  sizes  at  any  instant,  a  di.stribution  of  rise 
velocities  exists.  However,  for  simplicity  of  this  First  generation  model,  we  make  the  assumption  that 
the  number  of  crystals  is  a  constant  and  consider  the  mean  radius  /•  of  the  frazil  as  a  function  of  time: 


I  Af(r)  a 
bnNH  Pi) 


and  the  frazil  crystal  thickness  w  as  a  function  of  time: 


h-(t)  = 


2^ 

a 


where  a  =  a  representative  diameter/thickness  ratio  of  the  crystals 
Pi  =  the  density  of  ice 

N  =  the  number  of  crystals  per  cubic  meter,  typically  between  1  O  Vm-^  and  1 0^/m  ’ 


(4) 


(5) 
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We  assume  that  during  the  time  between  renewals  the  crystals  rise  at  their  terminal  velocity;  the 
velocity  at  which  the  buoyancy  force  balances  the  drag  force.  Wuebben  (1984)  observed  that  frazil 
crystals  in  a  quiescent  water  column  rise  without  rocking  or  tumbling  in  a  favored  orientation  of  basal 
plane  normal  to  the  vertical.  Based  on  this  force  balance,  and  assuming  the  plane  of  the  crystal  disk 
to  be  horizontal,  we  can  express  the  terminal  velocity  ^(t)  as 


V{x)  = 


4,?  V(T) 
.  aCd 


1/2 


where  g'  =  ,?(l-Pi/p) 

g  -  acceleration  due  to  gravity 
Cd  =  the  drag  coefficient. 


(6) 


Gosink  and  Osterkamp  (1983)  obtained  Cd  by  fitting  observed  data  as 
logCd=  1.386-0.8921  log  /? +0.111  (log/?P  (7) 


where  the  Reynolds  number  R  is  based  on  the  rise  velocity  and  the  diameter  of  the  frazil  disk.  They 
found  that  the  rise  velocity  calculated  from  eq  6  and  7,  with  diameter-to-thickness  ratios  of  1 0  and  50, 
bracketed  the  observed  values  for  disks  with  diameters  between  1  and  3  mm.  These  data,  presented 
in  Figure  1 ,  included  both  quiescent  laboratory  water  columns  and  agitated  water  columns  in  the  field. 
All  of  these  data  were  comparable.  The  scatter  may  be  a 
result  of  the  small  eddies. 


THE  MODEL 

Our  model  is  based  on  several  hypotheses.  The  scale  of 
the  energy-containing  eddies  is  comparable  to  the  depth  of 
the  shear  flow.  During  a  surface  renewal  the  frazil  is  mixed 
downward,  primarily  by  advection  by  the  large  eddies.  We 
neglect  the  mixing  caused  by  smaller  eddies.  Between 
surface  renewals  the  buoyancy  of  the  growing  crystals 
dictate  upward  movement.  The  frazil  crystals  at  the  sur¬ 
face  bond  if  the  concentration  is  high  and  the  time  of 
exposure  to  subfreezing  air  is  sufficiently  long.  The  for¬ 
mation  of  crystal  floes  at  the  surface  results  in  increased 
buoyant  force  and  resistance  to  reentrainment  by  the  large 
eddies. 

The  water  surface  is  occasionally  renewed  at  individual 
locations  by  the  large  eddies.  We  will  not  characterize  the 
mixing  action  of  these  eddies  individually.  Instead  we 
consider  their  gross  effect  on  the  frazil  concentration 
distribution  by  averaging  local  unit  areas  over  the  entire 
surface  with  a  probability  distribution  function  (PDF).  We 
propose  the  following  PDF  to  perfonn  the  averaging; 


Figure  1 .  Rise  velocity  of  frazil  crystals 
(afterGosinkandOsterkamp  1983 ).  The 
upper  and  lower  tines  represent  diameter! 
thickness  ratios  of  10  and  50,  respectively. 


F(t)  =  Je 


-ft 


1-e-^" 


(8) 


4 


Crystal  Layer 


Concentration 


Figure  2 .  Dimensionless  average  surface  age 
as  a  function  of  dimensionless  time. 


Figure  3.  Local  frazil  concentration  pro¬ 
file. 


F(0  is  positive  and  for  0  <  /  <  x  has  an  integral  of  1 ,  requirements  fora  PDF.  Dank  wert  s  (1951)  proposed 
a  similar  PDF  for  modeling  gas  absorption  across  a  free  surface,  and  eq  8  approaches  Dank  werts  ’  PDF 
as  T  becomes  large.  Komori  (1982)  experimentally  demonstrated  the  utility  of  Dankwerts’  PDF  by 
establishing  a  mass  transfer  coefficient  for  the  liquid.  We  chose  the  PDF  given  in  eq  8  because  the 
mixing  of  a  surface  layer  into  the  interior  of  a  flow  is  the  same  for  both  gas  absorption  and  frazil 
development. 

The  residence  time  of  frazil  at  the  surface  is  important  in  the  development  of  structure,  and  can  be 
characterized  by  surface  age.  The  age  of  a  surface  area  is  the  time  t  since  the  last  renewal  or  since  t^, 
whichever  is  shorter.  At  time  x  the  age  of  the  water  surface  patches  vary  between  0  and  x.  and  the 
fraction  with  ages  between  t  and  t  +  dtis  F(t)dt.  The  average  surface  age  7  is 

t'  »  f  <F(r) 


,9, 

1-e 

The  dimensionless  average  surface  age/7  as  a  function  of  dimensionless  time  fi  obtained  from  eq  9 
is  presented  in  Figure  2. 

If  the  rise  velocity  of  frazil  crystals  in  turbulent  water  is  known,  the  time  for  a  crystal  at  the  bottom 
of  the  shear  flow  to  rise  to  the  water  surface  t*  is 

t*=^i—.  (10) 


For  values  of  r  less  than  t*,  the  crystals  move  upward,  leaving  a  space  cleared  of  frazil  near  the  bottom. 
Meanwhile,  cr^'stals  accumulate  in  a  surface  layer  of  concentrated  frazil.  The  local  frazil  concentration 
profile  below  the  surface  layer,  shown  in  Figure  3,  follows  from  the  assumptions  of  uniform  crystal 
size  with  a  single  rise  velocity,  and  a  uniform  concentration  after  large  eddy  mixing.  This  profile  is 
expressed  as 


c(yd)  = 


_  (^)  for 


L>± 

H  t* 
-L 

H  t* 
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where  C  (x)  is  the  mean  volumetric  concentration  of  frazil  in  the  water  column.  C  (t)  depends  on  the 
rate  of  heat  loss  at  the  surface  Q  and  time  x.  For  t  >  t*  all  crystals  are  gathered  in  the  surface  layer. 

Frazil  crystals  that  gather  at  the  water  surface  and  are  undisturbed  for  a  sufficient  period  of  time  may 
become  bonded  and  resist  reentrainment  into  the  flow  interior  by  turbulent  eddies.  The  mean 
volumetric  concentration  of  frazil  in  the  water  column  after  can  be  expressed  as 

C(x)=p(x)!^  (12) 

Pi// 

where  |i(x)  represents  the  mass  fraction  of  frazil  in  circulation  relative  to  the  total  mass  of  frazil. 
Initially  all  the  frazil  crystals  are  in  circulation  and  P(0)  =  1.  Subsequently  P(x)  may  decrease  if 
structured  frazil  forms  at  the  surface  and  resists  reentrainment. 

Using  F(t)dt  as  the  weighting  factor,  we  obtain  an  average  or  global  vertical  concentration  profile 
c(y,x)  from  the  individual  profiles  under  areas  of  different  age  as 

c(y,  x)=  [  C{y,t)F(t)dt  (13) 


With  the  average  concentration  profile  we  can  establish  the  average  thickness  of  the  corresponding 
surface  frazil  layer  fi(x).  A  balance  equation  written  for  the  total  ice  volume  on  a  unit  surface  area  and 
in  the  water  column  below  it  at  time  x  yields 


M(x)-p,\  c(y,x)dy 


where  Q  is  the  frazil  concentration  of  the  surface  layer.  Martin  and  Kauffman  (1981)  obtained  the 
range  of  Cj  as  0.2-0.4.  Our  field  observations  of  frazil  in  rivers  indicate  a  lower  Q  value.  Data  from 
Tsang  (1988)  indicate  a  Cj  of  0.03  in  a  river.  The  maximum  thickness  of  the  surface  layer  can  be 
obtained  from  eq  14  as 

«  -W(t) 


We  visualize  a  critical  surface  layer  thickness  as  a  threshold  for  the  ice  at  the  surface  to  develop 
structure.  This  critical  thickness  depends  on  the  freeboard  of  crystals  at  the  surface  and  other 
environmental  parameters.  When  the  layer  thickness  exceeds  the  ice  protrudes  above  the  water 
surface  and  has  remained  there  for  a  sufficient  time  to  promote  the  bonding  of  crystals.  These  ice  floes 
have  increased  resistance  to  reentrainment  and  P(x)  <  1.  Before  this  thickness  is  attained,  all  crystals 
on  the  surface  are  entrained  back  into  the  interior  when  the  surface  is  renewed  and  P(x)  =  1 .  We  assume 
that  P(x)  is  a  function  of  the  average  surface  layer  thickness  and  the  average  surface  age,  and  propose 
a  simple  expression  that  is  correct  at  both  extremes  of  time  as 


p(x)  =  1  -  — 


fi(X)  = 


0  for 

B(x)-Bc  for  B(x)>Bc 


The  development  of  a  surface  layer  will  reduce  the  rate  of  heat  loss  Q(x)  from  the  turbulent  water 
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column  and  inhibit  the  growth  of  the  frazil  crystals  in  circulation.  When  Q(t)  -  0,  crystal  growth  stops, 
and  the  rise  velocity  ceases  to  increase. 


UNSTRUCTURED  FRAZIL 


We  can  examine  some  aspects  of  the  conceptual  mtxlel  using  the  idealized  local  concentration 
profile  shown  in  Figure  3.  Substituting  eq  8  and  1 1  into  13  and  performing  the  integration,  we  obtain 

~  for  T  >  ^  or  T//  *  > 

e(v.t)  ^  1  -  e/'^  '  H 

C(t)  _  . 

1  >  for  X  <  ^  or  x//  *  <  ^ 

H 

where  y/U  is  the  time  needed  for  a  crystal  to  ri.se  from  the  bottom  toy.  Given  a  specifc  depth  and  time 
the  relative  concentration  c(y,x)/C  (x)  depends  on  the  ratio  of  the  renewal  frequency  to  the  rise  velocity . 
We  use  a  dimensionless  renewal  frequency  ft*  to  characterize  the  relative  magnitude  of  the  downward 
motion  of  frazil  by  the  renewal  eddies  with  that  of  the  upward  motion  by  buoyancy  of  the  crystals. 
Figure  4  shows  the  mean  concentration  profiles  c(y.x)/ 

C(x)  through  dimensionless  depth  at  dimensionless  re-  too  io  i/t*  =  i 

newal  frequencies  of  0. 1 , 1  and  1 0,  and  several  values  of  VP  ’  ' 

dimensionless  time  x/r*.  We  note  in  Figure4  that  at  small  j-  ^^1/  - 

dimensionless  time  the  mean  concentration  is  C(x)  every-  ~05 

where  in  the  profile  except  near  the  bottom  of  the  shear  &  o.4  -jl  ^ — ■  - 

flow.  At  x/r*=  1  crystals  originally  at  the  bottom  arrive  at  0.2  ■jl  ^  — — 0.2 

the  surface.  If  a  local  renewal  has  not  occurred,  all  the  0  ’ 

crystals  reside  at  the  surface.  The  average  concentration  ^ 

profiles  at  x/r*  =  1  range  from  zero  at  the  bottom  to  C(x)  .  .  „ 

T/t  =  100. 10  1 

immediately  below  the  surface  layer.  ’  0 - 1 - ^ - 7 - 1 — 

At  larger  dimensionless  times  and  after  the  surface  o  s  -  tt*=  1  /  - 

heat  loss  has  ceased,  the  rate  of  frazil  entrainment  by  the  0  6  -  z' 

large  eddies  is  balanced  by  buoyancy  and  a  dynamic  g-  0  4  - 

equilibrium  develops.  A  portion  of  the  frazil  remains  in  *^02-  y  02 

suspension  indefinitely  beneath  a  surface  layer  of  con-  ^  - ; — ^  \  0,1 

stant  thickness  that  is  without  structure.  The  equilibrium  0  0.2  0.4  0.6  o.s  1.0 

frazil  concentration  profile  at  large  x  is  obtained  from  eq  Relative  Concentration 

17;  lOi - ^ ^ - 1 - 1 - i 


0.2  0.4  0.6  0.8 

Relative  Concentration 


^  ( 1  _  e-/-v  /'■)-(!-  e-A' *.'■  I")  (18) 

C(x) 

The  equilibrium  profile  with  dimensionless  depth  is  only 
a  function  of//*,  and  several  profiles  are  given  in  Figure 
5.  When  the  dimensionless  renewal  frequency  is  large, 
the  equilibrium  profile  reflects  well-mixed  conditions; 
when//*  is  small,  the  frazil  collects  in  a  layer  at  the  water 
surface.  Note  that  if  the  surface  heat  loss  persi.st,s  longer, 
the  mean  cry.stal  size  increases  and  ft*  becomes  smaller. 
Then,  with  all  other  conditions  the  same,  the  resulting 


vt*  =  100, 10. 1  - 


3.2  0.4  0.6  0.8 

Relative  Concentration 


Figure  4.  Relative  frazil  concentration 
c(y.T)/C(  T)  profiles  through  time  T/t’  for 
selected  dimensionless  surface  renewal 
freiptencies  ft*. 
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Figure  5  Equilibrium  relati ve  frazil  con¬ 
centration  c(yllC(T)  profiles  as  a  func¬ 
tion  of  dimensionless  surface  renewal 
frequency  ft*. 


Figure  6.  Nondimensional  equilibrium  surface  layer  thickness  as  a  function  of 
dimensionless  surface  renewal frequency.  The  curve  defines  regions  of  structured  and 
unstructured  frazil.  The  points  shown  correspond  to  the  cases  presented  in  Figures  8  and  9. 


equilibrium  state  would  have  generally  lower  relative  concentrations  throughout  the  depth  and  much 
more  frazil  in  the  surface  layer. 

An  equilibrium  surface  layer  corresponds  to  each  equilibrium  profile.  Substituting  eq  1 8  into  eq  1 4. 
performing  the  integration,  and  rewriting  using  eq  12  with  |3  =  1  and  eq  15  yields  the  equilibrium 
thickness  of  the  surface  layer  5^^  as 


^eq  _  ( 1  —  e  \ 
®  max  'ft  ' 


(19) 


The  dimensionless  equilibrium  surface  layer  thickness  obtained  from  eq  19  is  given  in  Figure  6. 
Equilibrium  thickness  approaches  fl^ax  ttsft*  approaches  zero,  and  it  approaches  zero  as  ft*  becomes 
large.  Thus,  without  bonding  between  the  crystals,  the  surface  layer  remains  unstructured,  with  an 
average  thickness  that  is  inversely  related  to ft*.  When ft*  is  sufficiently  small,  local  patches  of  crystals 
at  the  surface  become  thick  and  have  a  long  surface  residence  time.  The.se  conditions  promote  bonding 
between  crystals,  and  subsequent  renewal  eddies  may  not  entrain  all  the  frazil  on  the  patch  of  surface. 
The  partial  entrainment  that  results  is  accounted  for  in  the  model  by  P(t)  <  1 .  Bonding  initiates  the 
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formation  of  ice  floes  on  the  surface  and  a  continual  loss  of  crystals  to  the  surface  layer  that  will 
eventually  deplete  the  flow  interior. 


DEVELOPMENT  OF  STRUCTURED  FRAZIL 

Structured  frazil  begins  to  develop  at  the  water  surface  when  fl(T)  exceeds  the  thre.shold  represented 
by  Bf..  Beyond  this  threshold  the  equilibrium  thickness  no  longer  limits  the  layer  thickness,  which  will 
approach  Bmax  "''^h  time.  The  critical  surface  layer  thickness  made  nondimensional  by 
compared  with  the  dimensionless  equilibrium  thickness  for  the  same//*.  To  allow  the  possibility  of 
structured  frazil,  the  equilibrium  surface  layer  thickness  must  exceed  B,..  Therefore,  the  relative  values 
of  equilibrium  and  critical  layer  thickness  define  the  regions  of  structured  and  unstructured  frazil 
shown  in  Figure  6. 

We  obtain  a  set  of  equations  for  the  dimensionless  surface  layer  thickness  using  eq  8,  11. 

12,  13,  14  and  17  as 


fi(T)  _  1  -(l+/-)e-/^ 
fimax  //*(l-e-/") 


for  i  <  I  and  B  <  Be  , 
t* 


B(t)  _  I  -(e-f'  +ft*  e-f^) 

®max 


for  —  >  1  and  B  <  Br  , 


®  max 


l-e-/^(l+/'c) 

/  gc^gmax  j'j 

\1  —  gc^gmax 

[l_/r-e-/t(l+A -/'*)] 

/t*{l-e/^)| 

1  y/ 

1  [l_/,*_e-/x(l +/!-//*)] 

il-gc^gmax 

for  <  1  and  B>  Be ,  and 
t* 

1  -  {e-ff'  +  ft*  e-/^)  +  (  ^C^^max  y/j  [j  g-//*] 

g(^)  _  _ -  fic^gmax  / _ 

gmax  /f*(l  ~  e'^'')  +  I - L - ft\  [l -//*- e“^^*] 

\l-ec/gn,ax  I 


(20) 


for  ^  >  1  and  B>  Be, 
t* 

The  dimensionless  surface  layer  thickness  as  a  function  of  dimensionless  time  is  given  in  Figure  7  and 
the  mass  fraction  of  frazil  in  circulation  P  as  a  function  of  dimensionless  time  is  given  in  Figure  8  for 
ranges  of//*  and  B/B^gx-g  l  and//*>  10  the  surface  layer  remains  unstructured  with 

a  dynamic  equilibrium  between  the  frazil  on  the  surface  and  that  in  the  interior.  For  these  conditions 
P  is  constant  with  a  value  of  1 .  The  surface  layer  thickness  eventually  exceeds  B^  if  ft*  becomes  less 
than  10,  and  in  time  the  flow  interior  is  depleted  of  crystals.  B/B^gx  approaching  1  and  P  approaching 
zero  indicate  this  condition.  LargerB/6n,gx  values  allow  a  mixed  flow  dynamic  equilibrium  at  smaller 
dimensionless  renewal  frequencies.  The  final  surface  layer  thickness  can  only  be  Bn,gx  with  P  -  0,  or 
Bgq  with  P  =  1 . 
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Figure  7.  Dimensionless  surface  layer 
thickness  as  a  function  of  dimensionless 
time  for  selected  dimensionless  critical 
surface  layer  thicknesses  and  suiface 
renewal  frequencies. 


Fi}>ure  8.  pas  a function  of  dimension¬ 
less  time  Jor  selected  dimensionless 
critical  layer  thicknesses  and  surface 
renewal  frequencies. 


MODEL  APPLICATION 

We  now  compare  the  model  w'ith  empirical  relationships  and  field  observations  that  have  been 
reported  in  the  literature.  Matousek  (1984)  propo.sed  a  formula  for  the  critical  velocity  U  (in  m/s) 
required  to  maintain  well-mixed  frazil  in  a  river: 

t/  >  lii- iCfO.VC +6)]'^-  (21) 

1/2  ^ 

where  C  =  is  the  Chezy  coefficient  and  C.  is  the  dimensionless  Chezy  coefficient.  The  critical 
velocity  is  related  to  the  shear  velocity  as  U  -  u*C*.  Using  eq  21./=  3.73  tu/H,  and  a  constant  )  'of 
0.01  m/s  with  C*  values  of  5  and  30  yields  dimensionless  renewal  frequencies  of  19  and  16. 
respectively.  The  empirical  relationship  of  Gosink  and  Osterkamp  (1983)  for  approximating  the 
velocity  threshold  between  layered  and  well-mixed  frazil  isU-  C/42,  or  in  temis  of  shear  velocity  u. 
=  0.07.5  m/s.  Again,  assuming  \'  =  O.Ol  m/s  the  corresponding//*  =  28.  The  equilibrium  profiles  in 


lO 


1.2 


Figure  9.  P  and  B/B,^  as  functions  of 
normalized  dimensionless  time  for  se¬ 
lected  values  of  ft*  and 


B  and  D 


ft/ft 


max 


Figure  5  corresponding  to  these  empirical  criteria  all  represent  well-mixed  frazil.  In  the  model, 
equilibrium  profiles  vary  continuously  with  dimensionless  renewal  frequency.  However,  for  P(t)  <  1 
the  equilibrium  profiles  are  not  sustainable  and  the  flow  interior  is  progressively  depleted  of  crystals. 

Gosink  and  Osterkamp  ( 1983)  reported  field  data  on  the  Chatanika  River  and  Coldstream  Creek  in 
Alaska.  They  observed  well-mixed  frazil  in  the  former  and  layered  frazil  in  the  latter.  Using  the  mean 
depth  and  the  shear  velocity  values  reported,  and  again  taking  constant  V = 0.0 1  m/s.  the  dimensionless 
renewal  frequency  is  33  for  Chatanika  and  25  for  Coldstream.  Both  values  represent  modeled 
equilibrium  concentration  profiles  that  are  well  mixed.  The  discrepancy  between  these  results  and  the 
Coldstream  Creek  observations  may  be  caused  by  ice  on  the  water  surface  that  was  not  drifting  freely 
and  a  near-surface  velocity  that  was  much  slower  than  the  maximum  of  the  profile.  Since  the  equation 
for  renewal  frequency  assumes  a  logarithmic  velocity  distribution,  the  estimate  of ft*  may  not  be  valid. 

Based  on  the  renewal  frequency  relationship  and  the  assumed  constant  rise  velocity,  the  empirical 
relation.ships  and  field  observations  discussed  above  indicate  ft*  =  20  as  representative  of  the  flow 
conditions  at  the  threshold  between  mixed  and  layered  frazil.  The  model  provides  no  means  to 
determine  the  appropriate  value  of  6,.  or  which  must  come  from  laboratory  and  field 

measurements.  We  are  not  aware  of  available  measurements  that  would  provide  this  information. 
Using  Figure  6  we  obtain  BJB^^  «  0.05  for  the  model  to  reproduce  the  observations  of  mixed  or 
layered  frazil.  This  result  indicates  that  only  a  few  percent  of  the  total  mass  of  frazil  is  needed  at  the 
water  surface  to  initiate  bonding  between  crystals. 

We  will  now  apply  the  model  to  the  Chatanika  River  case  discussed  above  with/r*  and  BJB^^^  taken 
as  constants.  The  initial  mass  of  frazil  is  assumed  to  be  negligibly  small,  and  the  heat  loss  from  the  river 
is  assumed  to  be  constant  at  10  W/m^.  The  development  of  structure  in  the  surface  layer  depends  on 
^c/^max-  With//*  =  33,  layered  frazil  develops  iovBJB^^  =  0.03,  and  the  frazil  remains  well  mixed 
for  BJB^^  =  0. 1 0.  These  results,  expressed  as  P  and  B/Smax  time,  are  presented  in  Figure  9.  The 
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Figure  JO.  Normalized  frazil  concentration  as  a  function  of  dimensionless  depth  and  normalized  dimen¬ 
sionless  time  for  selected  values  of  ft*  and  Note  the  differences  in  the  concentration  scales. 


dimensionless  time  in  this  figure  ft  is  normalized  by  ft,^^  of  20.  The  equilibrium  surface  layer 
thickness  is  only  a  few  percent  of  Next,  we  reduce  the  shear  velocity .  corresponding  to  a  decrease 

in  the  intensity  of  the  turbulence,  and  obtain  ft*  =  5.  Layered  frazil  now  develops  at  =  0. 19. 

and  the  frazil  is  well  mixed  with  BJB^^^  =  0.30  (Fig.  9).  The  dimensional  time  represented  in  Figure 
9  is  much  greater  at  a  given  normalized  dimensionless  time  for  the  smaller  value  of ft*.  In  both  cases 
that  developed  layered  frazil,  the  surface  layer  approached  equilibrium  until  disrupted  with  P  <  I . 

A  set  of  three-dimensional  plots  of  normalized  frazil  concentration  in  depth  and  time  are  presented 
in  Figure  10  for  these  same  cases.  The  parameter  Cftmax^  P-0  used  to  normalize  the  concentration 
is  the  mean  volumetric  concentration  of  frazil  in  the  watercolumn  at  T,^ax  ‘f  assume  that  P  =  1 .  For 
ft*  =  33  the  concentration  profiles  at  early  times  are  vertically  well-mixed  and  normalized  concentra¬ 
tions  gradually  increase.  Subsequent!;,  the  collapse  or  continued  growth  of  the  profile  depends  on 
BJB^^.  At  ft*  =  5  the  normalized  concentration  is  not  uniform  in  the  vertical  direction,  but  all  other 
trends  are  the  same. 
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CONCLUSIONS 


A  model  is  presented  for  the  evolution  of  frazil  with  time  and  over  depth  in  a  turbulent  flow.  The 
mechanism  in  the  model  for  counteracting  the  buoyancy  of  the  crystals  is  that  large  energy-containing 
eddies  well  up  to  the  surface  and  entrain  all  or  a  portion  of  the  frazi  1  at  the  surface  back  into  the  interior. 
The  crystals  in  circulation  are  assumed  to  be  uniformly  distributed  over  the  depth  after  large  eddy 
mixing.  The  crystals  then  resume  upward  migration  at  their  turbulent  rise  velocity  until  interrupted  by 
the  next  renewal  eddy.  An  assumed  probability  density  function  is  used  to  obtain  averages  for 
parameters  over  an  ensemble  of  discrete  local  volumes  including  concentration  profile,  surface  age. 
and  surface  layer  thickness. 

Starting  from  a  well-mixed  initial  state  with  seed-size  crystals,  a  surface  layerof  frazil  develops  with 
time.  The  thickness  of  this  layer  depends  on  the  time  and  rate  of  heat  loss  through  the  water  surface, 
the  frazil  concentration  of  the  surface  layer  that  naturally  evolves,  and  the  dimensionless  renewal 
frequency.  This  frequency  is  the  product  of  the  frequency  of  surface  renewal  and  the  gathering  period 
of  frazil  to  the  surface.  The  gathering  period  depends  on  the  rise  velocity,  which  increases  with  time 
according  to  the  rate  of  heat  loss.  The  dimensionless  renewal  frequency  also  characterizes  the  frazil 
distribution  at  equilibrium  as  well-mixed  for  ft*  »  1 ,  or  buoyancy  dominated  and  layered  for  ft*<  1 . 
The  dimensionless  renewal  frequency  and  a  critical  surface  layer  thickness  determine  whether  the 
surface  layer  develops  structure.  Depending  on  the  development  of  structure  in  the  surface  layer,  the 
frazil  will  tend  toward  either  an  equilibrium  profile  and  surface  layer  thickness,  ora  zero  concentration 
profile  with  all  the  frazil  layered  at  the  surface.  Existing  empirical  criteria  and  field  data  on  the 
transition  between  layered  and  well-mixed  frazil  in  rivers  provide  a  consistent  surface  renewal 
frequency  at  the  transition  and  an  initial  estimate  of  the  critical  surface  layer  thickness. 

This  model  for  the  vertical  distribution  of  frazil  in  a  watercolumn  contains  parameters  that  quantify 
several  critical  processes  in  the  evolution  of  frazil  from  crystals  to  larger  ice  forms.  We  have  used  the 
model  to  systematically  assess  observations  of  mixed  and  layered  frazil,  and  more  detailed  analyses 
could  be  performed.  At  present  the  model  is  a  useful  tool  that  can  be  used  to  guide  future  observations 
of  frazil  evolution.  Additional  measurements  that  characterize  the  surface  renewal  frequency,  the 
turbulent  rise  velocity  of  frazil,  and  thecritical  surface  layer  thickness  are  needed  to  refine  and  provide 
further  validation  of  this  model. 
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